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ABSTRACT 

We present a joint gravitational lensing and stellar-dynamical analysis of 11 early-type galaxies 
(median deflector redshift Zd = 0.5) from Strong Lenses in the Legacy Survey (SL2S). Using newly 
measured redshifts and stellar velocity dispersions from Keck spectroscopy with lens models from 
Paper I, we derive the total mass density slope inside the Einstein radius for each of the 11 lenses. 
The average total density slope is found to be (7'} = 2.161q'q| (ptot °^ r ~ 7 )> with an intrinsic scatter 
of 0.251q'q7. We also determine the dark matter fraction for each lens within half the effective radius, 
i? e ff/2 and find the average projected dark matter mass fraction to be 0.42lg JJg with a scatter of 



0.20_QQ7 for a Salpeter IMF. By combining the SL2S results with those from the Sloan Lens ACS 
Survey (median z& — 0.2) and the Lenses Structure and Dynamics survey (median 2d = 0.8), we 
investigate cosmic evolution of 7' and find a mild trend d(^')/dzd = — 0.251q'i2- This suggests that 
the total density profile of massive galaxies has become slightly steeper over cosmic time. If this result 
is confirmed by larger samples, it would indicate that dissipative processes played some role in the 
growth of massive galaxies since z ~ 1. 

Subject headings: galaxies: fundamental parameters — gravitational lensing — 



1. INTRODUCTION 

Early- type {i.e. elliptical and lenticular) galaxies in the 
local universe are considered to be simple objects (e.g. 
Bertin & Stiavelli 1993; Merritt 1999; Ciotti 2009). In 
the central few kpc most of their mass is dominated by 
stars, while at larger radii there is convincing evidence - 
at least for the most massive systems - that dark mat- 
ter halos are dominant. Their stellar populations are 
relatively simple, dominated by old stars with little or 
negligible star formation (e.g. Renzini 2006). Remark- 
ably, many global properties, ranging from the chemical 
composition of their stars to their size, luminosity and 
mass of the central black hole, correlate tightly with their 
stellar velocity dispersion, a (e.g. Bernardi et al. 2005; 
Graves et al. 2009). 

In spite of this suspected simplicity, their formation 
and evolution are still poorly understood, and therefore 
they are the subject of many observational and theo- 
retical investigations. A number of observational facts 
have proved difficult to explain by theoretical models. 
These include: i) the tightness of the empirical correla- 
tions with a (e.g. , Bernardi et al. 2005; Graves et al. 
2009; Nipoti et al. 2009a); ii) the so-called downsizing 
trend of their stellar populations, i.e. the correlation be- 
tween mean stellar age and present day stellar mass (e.g. , 
Thomas et al. 2005; Treu et al. 2005; van der Wei et al. 
2005; Juneau et al. 2005; di Serego Alighieri et al. 2005); 
iii) the evolution of the upper end of their mass function 
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since z ~ 1 (e.g. , Bundy et al. 2005, 2007; Cimatti et al. 
2006; van der Wei et al. 2009; Hopkins et al. 2010); iv) 
the unusually compact size of high redshift massive red 
galaxies (e.g. , Treu et al. 1998; Daddi et al. 2005; van 
Dokkum et al. 2008; Saracco et al. 2009; Cassata et al. 
2010; Mancini et al. 2010; Newman et al. 2010). From a 
theoretical standpoint, it is clear that understanding the 
interplay between baryons, black holes, and dark matter 
is essential to develop a scenario that can quantitatively 
reproduce all observations. The physical processes that 
need to be accurately modeled in a successful theory ap- 
pear to include black hole accretion and the related en- 
ergy and momentum feedback (e.g. , Croton et al. 2006; 
Ciotti et al. 2009), dry and wet major mergers (e.g. , 
Khochfar & Silk 2006; Ciotti et al. 2007; Robertson et al. 
2006), as well as minor mergers and dry accretion of mi- 
nor satellites (e.g. , Naab et al. 2009; Hopkins et al. 2010). 

Most of the observational studies, including those 
listed in the previous paragraph, are concerned with 
global parameters of early-type galaxies. An entirely new 
line of investigation can be opened up if we are able to 
dissect early-type galaxies and map their internal dynam- 
ical structure as a function of cosmic time. By decompos- 
ing the internal mass distribution of early-type galaxies 
into luminous and dark components we can start address- 
ing the following questions, the answers to which would 
provide essential clues as to their formation and evolu- 
tion. How and when is mass assembled to form early- 
type galaxies? How are baryons converted into stars and 
accumulated inside dark matter halos? Is the mass den- 
sity profile of early-type galaxies comparable to that ob- 
served in numerical simulations? Do isolated early-type 
galaxies undergo internal structural and dynamical evo- 
lution? As numerical simulations of early-type galaxies 
become more and more realistic, detailed knowledge of 
their internal structure (e.g. their distribution functions) 
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will provide more and more stringent tests of the current 
paradigm of structure formation (e.g. Meza et al. 2003; 
Naab et al. 2007; Onorbe et al. 2007; Lackner & Ostriker 
2010). 

Great progress in answering these questions has been 
achieved in the past few years with the systematic 
study of early-type galaxies acting as strong gravitational 
lenses. For these systems, strong lensing provides an 
absolutely calibrated measurement of mass at a fiducial 
radius (the Einstein radius), which is typically compa- 
rable in size to the effective radius. By combining this 
mass tracer with traditional diagnostics such as stellar 
velocity dispersion (Miralda-Escude 1995; Natarajan & 
Kneib 1996; Treu & Koopmans 2002), and stellar mass 
maps from multicolor imaging and/or spectroscopy, one 
can break many of the degeneracies inherent to each 
method alone, including the mass-anisotropy degener- 
acy, bulge- halo degeneracy, and the stellar mass/initial 
mass function (IMF) degeneracy (e.g. Koopmans & Treu 
2003; Treu & Koopmans 2004; Treu et al. 2010). Ad- 
ditional information can be gathered with the addition 
of weak-lensing (Gavazzi et al. 2007; Jiang & Kochanek 
2007; Lagattuta et al. 2010), although at the moment 
this is not possible for individual galaxies. The SLACS 
team applied this methodology to a sample of more than 
80 early-type galaxies (Bolton et al. 2006, 2008a; Auger 
et al. 2009), which have been shown to be indistinguish- 
able from equally massive non-lensing early-type galax- 
ies in terms of their internal properties and environment 
(Bolton et al. 2006; Treu et al. 2006, 2009). 

Among the most relevant findings of SLACS is that 
the total mass-density profile p to t °c 01 early-type 
galaxies is close to isothermal with 7' = 2.085^°'^ with 
intrinsic scatter less than 0.1 (Koopmans et al. 2009; 
Barnabe et al. 2009; Auger et al. 2010), even though 
neither the stars nor the dark matter obey a simple 
power law profile (see also Wucknitz et al. 2004). This 
"bulge-halo conspiracy" , similar to the disk-halo conspir- 
acy found for spiral galaxies (e.g. van Albada & Sancisi 
1986), has implications both for lensing studies but also 
for galaxy formation studies. Since the total mass den- 
sity profile is preserved by dry-mergers (Dehnen 2005; 
Kazantzidis et al. 2006; Nipoti et al. 2009b), and dark 
matter-only profiles are not isothermal (e.g. Navarro 
et al. 2010), the isothermal nature has to be established 
through dissipational processes (e.g. , Koopmans et al. 
2006). The other main finding of the SLACS survey is 
that (assuming a constant stellar IMF) the fraction of 
dark matter /dm within a fixed fraction of the effective 
radius increases with galaxy mass or stellar velocity dis- 
persion (see also Jiang & Kochanek 2007; Grillo et al. 
2008; Cardone et al. 2009; Cardone & Tortora 2010). 
Possible explanations for this trend are: varying effi- 
ciency in converting baryons into stars as a function of 
halo mass, varying inner slope of the dark matter halo 
with mass, and evolutionary processes. For example, dry 
mergers can increase the fraction of dark matter within 
the effective radius (e.g. , Nipoti et al. 2009b). 

The SLACS, however, sample is limited to low red- 
shift by the selection function of the parent sample of 
SDSS luminous galaxies. Therefore, evolutionary stud- 
ies with the SLACS sample are limited to a short base- 
line. Suitable samples of strong lenses are smaller at 



high rcdshift. In fact, most of the strong lensing galaxies 
known to date at z > 0.4 arc too faint and/or dominated 
by strongly lensed bright quasars to allow for detailed 
kinematic studies of the deflector. Of the handful of ex- 
ceptions (e.g. , Faure et al. 2008; Lagattuta et al. 2010), 
only a few of them have published lensing and dynamical 
analysis (Treu & Koopmans 2002; Ohyama et al. 2002; 
Treu & Koopmans 2004; Suyu et al. 2010). These early 
studies found results similar to SLACS, i.e. the internal 
slope 7' is close to isothermal, and the central dark mat- 
ter fraction appears to increase with mass as found by 
the LSD Survey (Treu & Koopmans 2004). However, the 
results also hint that things may have been different at 
z ~ 1 when the universe was less than half its present 
age: the scatter in 7' might have been larger (Treu & 
Koopmans 2004) and perhaps even the average might 
have been different (Koopmans et al. 2006). Unfortu- 
nately, current samples beyond z ~ 0.4 with lensing and 
dynamical data are too small to probe the evolution of 
the internal structure of early-type galaxies, its scatter 
and trends with mass at the same time. 

In this paper we present the first detailed study of a 
sample of early-type lens galaxies identified by the Strong 
Lenses in the Legacy survey (SL2S). We use new spec- 
troscopic data obtained at the W.M. Keck telescopes, in 
combination with multicolor photometry from the CFHT 
Legacy Survey to infer dynamical and stellar mass for 
the deflector galaxies. This imaging data-set is comple- 
mented by EST imaging that allowed to confirm the lens- 
ing nature of the SL2S candidates and detailed lens mod- 
eling. The newly measured source and deflector redshifts 
are combined with gravitational lens models from Paper I 
(Gavazzi et al., in preparation) to infer lensing masses. 
We determine the total mass density profile of the early- 
type lens galaxies (quantified by 7') and their central 
dark matter fraction (/dm)- The median redshift of the 
sample is z = 0.494, providing an ideal complement to 
the earlier SLACS and LSD samples. By combining the 
three samples we investigate evolutionary trends in these 
quantities. 

This paper is organized as follows. In Section 2 we 
describe the SL2S survey, and how the galaxy-scale lens 
candidates were selected. We then present our spectro- 
scopic measurements in Section 3 before combining them 
with lens model parameters to investigate the mass struc- 
ture of massive galaxies since redshift 0.9 in Section 7. 
In Section 4 we incorporate stellar mass estimates in or- 
der to separate the dark and luminous components of 
the lens galaxies. The cosmic evolution of the total mass 
density slope and dark matter fraction are discussed in 
Section 8. After a brief discussion of our results in Sec- 
tion 9 we conclude in Section 10. Throughout this pa- 
per magnitudes are given in the AB system. We assume 
a concordance cosmology with matter and dark energy 
density Q m = 0.3, f^A = 0.7, and Hubble constant H o =70 
km s _1 Mpc _1 . 

2. STRONG LENSES IN THE LEGACY SURVEY: SL2S 

In this section we describe our feeder survey, the CFHT 
Legacy Survey, and how we identified our new sample 
of galaxy scale lenses at intermediate redshift. A more 
detailed description is given in Paper I. 

Strong lensing candidates were selected from the 
Canada-France-Hawaii Telescope Legacy Survey 
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(CFHTLS) 7 . In brief, the survey consists of two main 
components of sufficient depth and image quality to be 
interesting for lens searching 8 . Both are imaged in the 
u*, g, r, i and z bands with the 1 deg 2 field-of-view 
Megacam Camera. The multi epoch Deep survey covers 
4 pointing of 1 deg 2 each. Two different image stacks 
were produced: D-85 contains the 85% best seeing im- 
ages whereas the D-25 only includes the 25% best seeing 
images. For finding lenses we only considered these 
latter, better resolution stacks. They reach a typical 
depth of u* ~ 26.18, g = 25.96 ~ 25.47, r ~ 25.43, 
i ~ 25.08 and z ~ 24.57 (80% completeness for point 
sources) with typical FWHM point spread functions of 
0'.'75, Of! 69, 0'.'64, 0'.'62 and 0'/61, respectively. The Wide 
survey is a single epoch imaging survey, covering some 
171 deg 2 in 4 patches of the sky. It reaches a typical 
depth of u* ~ 25.35, g ~ 25.47, r ~ 24.83, i ~ 24.48 and 
z ~ 23.60 (AB mag of 80% completeness limit for point 
sources) with typical FWHM point spread functions of 
0'.'85, 0'.'79, 0'.'71, 0'/64 and 0'.'68, respectively. Because 
of the greater area, the Wide component is our main 
provider of lens candidates. 

Images from both the Deep and Wide survey modes 
were analyzed to find strong lens candidates using sev- 
eral algorithms, as described by Cabanac et al. (2007) at 
the group and cluster mass scales, and by Gavazzi et al. 
(in preparation) at the galaxy scale with the RingFinder 
algorithm. The ring-detecting algorithm is aimed at de- 
tecting compact rings around centers of isolated galaxies 
(< lO 13 ft- _1 M ), and works by focusing on the achro- 
matic image excesses around early-type lens galaxies that 
are indicative of the presence of lensed arcs. For each of a 
sample of pre-selected bright (zab < 22.5) red galaxies, a 
scaled, PSF-matched version of the i-band cutout image 
was subtracted from the g-band image of the same sys- 
tem. The rescaling in this operation is performed such 
that the early-type galaxy light is efficiently removed, 
leaving only objects with an SED different from that of 
the target galaxy. These (typically) blue residuals are 
then characterized with an object detector, and analyzed 
for their position, ellipticity, and orientation, and those 
showing characteristic properties of lensed arcs are kept 
as lens candidates. A sample of several hundred good 
candidates were visually inspected, and ranked for follow- 
up with HST. 

Currently, 65 CFHT galaxy-scale lens candidates have 
been observed firstly with ACS, then WFPC2, and fi- 
nally with WFC3 as snapshot programs over cycles 15, 
16 and the ongoing cycle 17. Details of these observa- 
tions are given in Paper I. Approximately 50% of the 
lens candidates were confirmed as lenses in this way. The 
sources are all faint blue galaxies, with very few showing 
signs of an active nucleus. Those with the most convinc- 
ing lens models (see Paper I) were selected for spectro- 
scopic follow up to obtain high precision redshifts for lens 
and source galaxies, and lens galaxy velocity dispersions. 
These observations are described in the next section. 

3. SPECTROSCOPIC OBSERVATIONS 

Follow up spectroscopy of 17 selected lens candidates 
was obtained using the Low Resolution Imaging Spectro- 

7 See http://www.cfht.hawaii.edu/Scicnce/CFHLS/ and links 
therein for a comprehensive description 

8 http : //terapix . iap . f r/rubrique . php?id_rubrique=259 



TABLE 1 

Observing logs 



run 


obsdate 


graname 


grisname 


slit 


time 


seeing 


1 


20 Jul 2006 


831/8200 


600/4000 


l'/5 


20 


tf'84 


2a 


23 Dec 2006 


831/8200 


300/5000 


1'.'5 


60 


0'.'60 


2b 


23 Dec 2006 


831/8200 


300/5000 


1'.'5 


40 


0'.'80 


3a 


13 Sep 2007 


400/8500 


400/3400 


T.'O 


45 


0'.'73 


3b 


13 Sep 2007 


400/8500 


400/3400 


T.'O 


30 


Of.' 7 5 


4a 


14 Sep 2007 


400/8500 


400/3400 


l'.'O 


30 


0"55 


4b 


14 Sep 2007 


600/7500 


600/4000 


Of.' 7 


30 


0"55 


4c 


14 Sep 2007 


831/8200 


600/4000 


l'.'O 


90 


0'.'b5 


5 


9 Sep 2009 


600/7500 


300/5000 


l'.'O 


30 


l'.'O 


6 


14 Jan 2010 


600/7500 


300/5000 


l'.'O 


45 


1'.'3 



Note. - The exposure time is given in minutes. The 
plate scale for both the red and blue chips is O'.'135/pixel for 
observations from 2009 and after. For observing runs before 
2009, the red chip had plate scale of 07211/pixel. 



graph (LRIS) on the Keck I telescope over six nights from 
2006 to 2010, with the aim of measuring the deflector and 
source redshifts and the velocity dispersion of the deflec- 
tor. Because of incompleteness in the execution of the 
HST snapshot observations mentioned above, we have 
spectra of 6 systems where -ffSTdata is not yet available. 
These lens candidates, however, were very promising on 
the basis of their CFHT images alone. In this section 
we describe the observations and their analysis in some 
detail, outlining our methodology for measuring the lens 
redshifts, source redshifts, and lens velocity dispersions. 

3.1. Details of the observations and data reduction 

Observations were made in long slit mode, with the slit 
centered on the deflector. The slit orientation was chosen 
to maximize both the source flux and the spatial separa- 
tion of the deflector and arc. The total flux is dominated 
by the deflector, so maximizing the source flux maximizes 
the likelihood of measuring the source redshift, z s . Spa- 
tial separation within the slit also separates the deflector 
and arc traces in the 2D spectra, allowing us to look for 
source emission lines directly, as discussed in Section 3.4. 
The seeing ranged between Of! 55 and l'/3 and total expo- 
sure times varied from 20 to 45 min, typically with two 
consecutive 15 min intervals. Between each exposure, we 
dithered along the slit by 10" to improve defect removal 
and sky subtraction. A summary of each observing run 
is given in Table 1. 

The later setups took greater advantage of the blue 
sensitivity of LRIS. We found that the most effective 
method to measure all three quantities was using a 680 
dichroic, so that the deflector redshift and velocity dis- 
persion could be measured in the blue. The 680 dichroic 
allowed us to choose a large wavelength coverage for the 
red, which was useful in detecting [On] out to z > 1.1. 

The data was reduced using a Python pipeline (de- 
veloped by MWA). On the red side, the night sky lines 
were used to determine the wavelength solution, while 
standard arclamps were used on the blue side. Two dif- 
ferent extraction windows were chosen: one wide extrac- 
tion window to ensure that light from the source was also 
included in the final spectrum (half widths between 5 and 
12 pixels from the central trace) and another narrower 
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SL2SJ140614+520252 
z d = 0.4797 

z = ? 



SL2SJ140123+555705 
z d = 0.5263 
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z a = 0.7044 
z = ? 
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6500 7000 
Wavelength (A) 



Fig. 1. — Keck/LRIS spectra for 17 SL2S lenses. An elliptical galaxy template is shown in a thin black line beneath each spectrum. The 
dotted black lines mark features of the lens galaxy, while source emission lines are highlighted in red. The shaded regions indicate A and 
B band telluric absorption features. 



extraction window to increase the signal to noise ratio, 
which is important for measuring velocity dispersions. 
The reduced spectra are shown in Figure 1. 

3.2. Measuring deflector redshifts 

Deflector redshifts were measured using the centroids 
and known rest frame wavelengths of prominent ab- 
sorption features. In the majority of cases, a mini- 
mum of the Call H and K lines were used, however 
other absorption features were also used when avail- 
able. The additional absorption features used were: 
Rn (3835A), G band (4305A), H 7 (434lA), H/3 (486lA), 
Mgll (5175A), Nan (5892A) and Ha (6563A). Typi- 
cally, between three and five absorption features were 
centroided to measure the redshift. The measured red- 
shifts are listed in Table 2. 

3.3. Measuring deflector velocity dispersions 

Velocity dispersions were measured by fitting combi- 
nations of stellar spectra over regions with prominent 
absorption features and high signal to noise. Linear com- 
binations of stellar spectra were used to fit a model to the 
data and calculate a velocity dispersion. A Python based 
implementation of the van der Marel (1994) velocity dis- 
persion code, developed by MWA and described by Suyu 
et al. (2009), was used. We use a set of templates from 



the INDO-US stellar library containing spectra for a set 
of seven K and G giants with a variety of temperatures 
and spectra. K and G giants were used because they pro- 
vide a good description of the spectra of our deflectors, 
as expected for massive ellipticals. 

The value of the velocity dispersion, a for each deflec- 
tor was determined by finding a consistent value over 
several spectral regions and features. If the mean signal 
to noise (S/N) ratio per rest- frame angstrom was ^ 10 
in the 4000-5000A range, then consistent and reliable a 
values could be measured. Of the 17 SL2S lenses, 12 
spectra had sufficient S/N to measure a. 

In general, the rest frame 4000-5000A range was used, 
as the G-band absorption feature is often uncontami- 
nated by atmospheric absorption, and contains no sharp 
changes in the continuum and the CCD efficiency is also 
good over this range. The average signal to noise ratio 
per angstrom in the 4000-5000A range for the 12 objects 
with measured velocity dispersions is 24.8. 

Other regions that were used to fit the stellar templates 
were: the 5000-6000A range (matching the Mgll and 
Na II absorption features) and the 3500-4000A range was 
also used, but only where the depth and continuum fit 
to the spectra were good. Generally, three regions of the 
spectrum that produced good fits and consistent a val- 
ues were used to calculate the final o and its associated 
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SL2SJ021737-051329 
z d = 0.6458 
z, = 1.847 



-+- 
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SL2SJ022648-040610 
z d = 0.7663 
z« = ? 



SL2SJ021411-040502 
z a = 0.6080 

z_ = ? 
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Wavelength (A) 
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Fig. 1. — continued. 



uncertainty. For each lens, one of the models generated 
to measure a is shown in Figure 2. Regions where atmo- 
spheric absorption was a problem, or the templates did 
not produce a good fit were masked out, as shown in the 
grey regions of Figure 2. The results of the velocity dis- 
persion measurements and the mean S/N per angstrom 
in the rest frame 4000-5000A range for each object are 
listed in Table 2. The mean measured a for the SL2S 
sample is 250 km s . The measured velocity disper- 
sions were then corrected to a uniform physical aperture 
using the slit width, the size of the extraction window 
and the empirical power-law relation of J0rgensen et al. 
(1995). The corrected velocity dispersion, er e 2, measures 
the velocity dispersion at i? e ff/2, as used by Bolton et al. 
(2008b). Inside the effective radius, the relation is well 
described by a power law: 

n \ -0.04 



and y are the 



2 7"ap , 

where 2r ap ss 2(xy /it) 1 / 2 and where 
width and length of the rectangular aperture. 

3.4. Measuring source redshifts 

Source redshifts were measured for five objects. Two 
objects, SL2SJ021737-051329 and SL2SJ141137+565119 
(z s = 1.847 and 1.420, respectively) have multiple source 
emission lines that were centroided to measure the red- 
shift. As can be seen in Figure 1, SL2SJ021737-051329 



has narrow emission lines typical of Type II AGN 
(Civ A1549, Hen A1640, Oni] A1666 and Cm] A1909). 
SL2SJ141137+565119 shows a clear splitting of the 
[O II] doublet as well as strong C in] and O in] emission. 

The three remaining source redshifts were identified 
using [On] AA3726.1, 3728.8 only. SL2SJ022610-042011 
and SL2SJ022511-045433 [z s = 1.232 and 1.1988, respec- 
tively) show a clear splitting of the doublet in the 2D 
spectra. For the objects where there was no clear split- 
ting of the [O n] doublet, the redshift measurement was 
more difficult. We used two methods to search for source 
emission lines. Firstly, we looked for emission lines in the 
2D spectra at the expected position of the source trace. 
Secondly, we looked at the residuals from both a simple 
fit to a standard elliptical galaxy template and also the 
residuals from the velocity dispersion measurement, dis- 
cussed in Section 3.3. Despite that three source redshifts 
were identified using [O n] alone, no other strong emission 
features are observed blue-ward of this feature, corrobo- 
rating the identification, as we would expect to see other 
lines if the detected feature were redder (e.g. Oni], Ha 
or H/3). 

The 2D spectrum of SL2SJ022511-045433 shows a 
slight offset between the [O n] emission of the arc and 
the arclet. We estimate the difference in redshifts to 
be Az = 0.0005, corresponding to a relative motion of 
~150 km/s. This is most likely due to kinematic struc- 
ture of the source. Consistent with this interpretation, 
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Fig. 1. 



continued. 



the lens model indicates the presence of two separate 
peaks in the surface brightness distribution of the lens, 
i. e. possibly a pair of source galaxies (see Paper I for de- 
tails). 

3.5. Systems with no source emission lines 

For the remaining 11 sources, for which no source emis- 
sion lines could be detected, we must estimate the source 
rcdshift from whatever information we have. We first 
consider the the HST and CFHTLS photometry. We 
did not attempt to infer photometric redshifts, since the 
disentanglement of lens and source colors in the low- 
resolution, multi-filter CFHTLS data was deemed likely 
to lead to significant systematic uncertainty. Instead, we 
conservatively used the redshift distribution of the faint 
galaxies in the COSMOS survey to provide a broad prob- 
ability density function (PDF) for z s . 



Leauthaud et al. (2007) model this distribution using 
the following functional form: 



Pr(z s |m) cx zl exp 



zo(m) 



1.5 



z (m) ■■ 



(0.18 m- 3.3) 
1.412 ' 



(2) 



(3) 



Here m is the AB magnitude of the source in the F814W 
filter. In most cases we have either a F606W magni- 
tude from the HST image modeling, or a g-band magni- 
tude from the CFHTLS image modeling: we assume that 
the sources, as faint blue galaxies, have spectral energy 
distributions consistent with flat in the AB magnitude 
system, and hence substitute our blue unlensed source 
magnitudes (output from the lens inversions) directly. 
This approximate transformation introduces a small ad- 
ditional uncertainty which we neglect relative to the in- 
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Fig. 1. — continued. 



trinsic PDF width. 

This model source redshift distribution was used by, 
among others, Gavazzi et al. (2007) and Lagattuta et al. 
(2010) when estimating the redshift distribution of back- 
ground sources. We note that this approach is qualita- 
tively different from that taken by Gavazzi et al. (2007) 
and Lagattuta et al. (2010): in their weak lensing stud- 
ies they did not have a single well-defined source mag- 
nitude, but instead integrated over the number counts 
down to the magnitude limit. Whereas we are able to 
use the small amount of information that the source 
brightness provides. We neglect the small uncertainty in 
zo(m) and the photometric uncertainty in m, and trun- 
cate the density to zero at z s < z&. We then draw sam- 
ples from Pr(z s |m) that we can then transform into dis- 
tances, physical masses and so on. This gives the broader 
PDF, Pr(z s ). 



While the source magnitude provides a very rough pho- 
tometric redshift, we can also use the lens geometry to 
give a similarly rough geometric redshift. This involves 
multiplying the COSMOS prior PDF Pr(z s |m) by an ad- 
ditional distribution describing our prior knowledge of 
the lens strength, as follows. In practice we do this by 
importance sampling the COSMOS prior (see e.g. Lewis 
& Bridle 2002; Suyu et al. 2010, for descriptions of how 
this process works). Given that the velocity dispersion 
of the dark matter is approximately equal to the central 
stellar velocity dispersion, and the apparently universal 
(approximate) isothermal profile of lens galaxies (see e.g. 
Kochanek 1994; Koopmans et al. 2006), we can down 
weight some predicted z s ,pdf values based on the veloc- 
ity dispersions they predict. The total deflector mass is 
modeled as a Singular Isothermal Ellipsoid (SIE), with a 
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velocity dispersion 

a sm = 186.2 kms- 1 x I lE^£^l (4) 
V arcsec D s 

where D^s and D s are the angular diameter distances be- 
tween the deflector and source and observer and source, 
respectively. The calculated ctsie for each z SiP df was used 
to give aPr(cr S iE)- A joint prior, Pr(z s |m)Pr(cr S iE), using 
the additional information from ctsie was then calculated 
and used to tighten the constraint on z s . 

To calculate Pr(crsiE), we assumed that the SL2S lens 
galaxies are at the high mass end of the velocity function 
of bright galaxies - high mass, because we know they 
are acting as gravitational lenses. To describe the veloc- 
ity function of the SL2S lenses, the velocity dispersion 
fitting function determined by Sheth et al. (2003) was 
used. The parameters in the fitting function were de- 
termined using measurements of a large sample of early- 
type galaxies drawn from the SDSS database. The fitting 
function was multiplied by cr 4 to mimic the lensing se- 
lection function (see e.g. Auger et al. 2010). However, 
the lenses, particularly those with measured c, are also 
luminosity selected because SDSS is a flux-selected sam- 
ple and more luminous objects are drawn from a large 
volume (Hyde & Bernardi 2009; Auger et al. 2010). To 
account for the luminosity selection, the Sheth et al. fit- 
ting function was heuristically multiplied by cr 4 . Fig- 
ure 4 shows the distribution of the measured SL2S ve- 



locity dispersions with the Sheth et al. fitting function, 
multiplied by two different selection functions, overlayed. 
The dashed curve shows the fitting function multiplied 
by <7 4 (mimicking the lensing selection function), while 
the solid curve shows the fitting function multiplied by 
both the lensing and luminosity selection functions. De- 
spite the small number of objects, the mean a of the 
SL2S deflectors (250 ± 39 km s _1 ) is in very good agree- 
ment with the peak of the solid curve (243 ± 48 km s _1 ). 
The prior adopted is then the Sheth et al. fitting func- 
tion multiplied by cr 8 ; the resulting distribution was then 
normalized and used to compute importances to weight 
the z s ,pdf values drawn from the COSMOS prior. We 
note that final result is relatively insensitive to the pre- 
cise choice of the exponent, since the dominant effect is 
the exponential cutoff, which eliminates unrealistically 
high stellar velocity dispersions. 

Our final prior PDF for the velocity dispersions of the 
SL2S galaxies is then given by: 

Pr(asiE) occr tt+7 e-(^)", (5) 
where the best fit parameters determined by Sheth et al. 
(2003) for early-type galaxies are: ct*=88.8±17.7 km s _1 , 
a=6.5±1.0 and /3=1.93±0.22. To account for the scat- 
ter in /gig = ctsie/c*, the function was conservatively 
smeared by 20%. This distribution is quite broad: the 
assumption of an approximately isothermal mass profile 
does not significantly bias our results when inferring the 
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Fig. 2. — Measuring the stellar velocity dispersions of the lens galaxies, by accurate absorption line fitting. For each object, we show 
the relevant portion of the LRIS spectrum (black line), compared to a model generated from all 9 INDO-US templates and a sixth order 
polynomial continuum (in red, with the continuum alone shown in green). The grey shaded wavelength ranges were not included in the 
fits. The lower sub-panels show the fit residuals in each case. 
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Fig. 3. — The [On] doublet is shown for each object with a measured source redshift and observed [On] doublet. 
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Fig. 4. — The distribution of measured velocity dispersions is 
shown with the Shcth et al. fitting function ovcrlayed. The solid 
and dashed curves show the fitting function multiplied by ct 8 (lens- 
ing and luminosity selection function) and a 4 (only lensing), re- 
spectively. The fitting functions were normalized to the unit area 
of the histogram. Note that the mean of the SL2S a distribution 
is in good agreement with the peak of the solid curve. 

profile slopes of individual lenses, since we include plau- 
sible scatter on a. Distributions of z s ,pdf for the final 
sample of lenses are shown in Figure 5. The dashed line 
shows z s ,pdf distributions from the COSMOS prior with- 
out weightings and the solid line indicates the final z S!p( jf 
for each of the lens. 

The primary effect of Pr(crsiE) is to disfavor lower red- 
shift solutions to the lens equation, that are offered by 
the COSMOS redshift distribution with its relatively low 
z s peak, which would indicate unrealistically large stel- 
lar velocity dispersions, well above 400 km s~ x . For each 
non-spectroscopic case, we give the median of each lens 
system's source redshift PDF in Table 2, along with the 
16 th and 84 th percentiles as indicators of the uncertainty. 
The method was verified by comparing the 5 spectro- 
scopically measured redshifts to those inferred via this 
procedure, using the same 5 systems. The grey verti- 
cal lines in Figure 5 show the measured source redshifts, 
which are in good agreement with z s .pdf after weighting 
by the Sheth et al. (2003) fitting function and the selec- 
tion function. 

4. PHOTOMETRIC MEASUREMENTS 

We now turn to our photometric data, in order to mea- 
sure sizes and stellar masses of our lens galaxies. The 
latter are required in order for us to probe the dark mat- 
ter fractions in the core regions of these objects, and are 
obtained by stellar population analysis of the CFHTLS 
photometry. 

4.1. Model size and magnitude estimation 

We fit an elliptically symmetric de Vaucouleurs profile 
(de Vaucouleurs 1948) surface brightness distributions to 
each of the lens galaxies in our sample, to measure the 
total magnitudes and effective radii. For each lens, we 
used the galf it software package to fit for the position, 
effective (half-light) radius, ellipticity, orientation angle 
and total magnitude. 



When available, we also used the HST data that gen- 
erally provide a better estimate of the deflector effec- 
tive radius. In particular, the light profile in the reddest 
available filter should be a better tracer of the underlying 
stellar mass distribution profile. Otherwise we used the 
average effective radius of the CFHT r, i and z bands and 
the internal filter-to-filter r.m.s. scatter to get a precise 
estimate of the effective radius and its associated mea- 
surement error. Under the assumption of a spheroidal 
distribution of stars, with a single stellar population and 
no color gradients, the scatter across the filter set in each 
morphological parameter gives a simple estimate of the 
uncertainty of each parameter, and in particular on i? e ff- 
In most cases, the typical error on R e g is of order 5%. 
Note also that i? e ff is expressed along the geometric mean 
axis, that is the radius of the circle enclosing the same 
area as the elliptical isophote enclosing half the light. 

To build up a picture of each lens' spectral energy dis- 
tribution (SED), with the deflector photometry devoid 
of light coming from the lensed background source, we 
follow a two-step iterative process to mask out the lensed 
images: first we fit models to the images with no mask- 
ing, and then use the residual image to manually mask 
out the lensed features, before refitting the model. Ap- 
parent magnitudes are then corrected for Galactic dust 
extinction (Schlegel et al. 1998). Errors are dominated 
by systematic photometric zero point uncertainties, dust 
correction errors and variations of galf it results de- 
pending on the masking strategy. All together these 
amount to systematic errors of ^0.05 magnitudes in each 
CFHT bands. The resulting effective radii and magni- 
tudes are given in Table 2. 

4.2. SED fitting methodology 

Estimates of the stellar masses were calculated using 
the CFHT galf it model magnitudes and composite stel- 
lar population models using a code developed by Auger 
et al. (2009). The code employs a Bayesian exploration of 
the stellar populations of galaxies using composite stel- 
lar population models produced by Bruzual & Chariot 
(2003). The code takes a set of photometric data for each 
object (g, r, i, z) magnitudes and their uncertainties from 
the CFHT filters and the measured redshift. The free pa- 
rameters of the star formation history are: the time when 
star formation began, the time scale of the exponential 
decay of the star formation rate, internal reddening due 
to dust extinction, and metallicity. The parameter space 
was then explored using a Markov Chain Monte Carlo 
(MCMC) routine, allowing a full determination of the 
posterior probability distribution function for each pa- 
rameter; see Auger et al. (2009) for a comprehensive de- 
scription of the method. We employ uniform priors for 
all of the model parameters (including the metallicity) 
due to the absence of any well-calibrated priors for the 
stellar population model parameters of massive galaxies 
at the redshifts of these lenses. 

The code was run for both Salpctcr (1955) and 
Chabrier (2003) initial mass functions (IMFs). The 
masses derived from the Salpeter IMF were a factor of 
~1.7 greater than the masses derived using the Chabrier 
IMF. In a recent study of 56 gravitational lenses iden- 
tified by the SLACS survey, Treu et al. (2010) found 
that for massive early-type galaxies, the Salpeter IMF 
provides stellar masses in approximate agreement with 
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Fig. 5.— The 

z s,pdi distributions are shown for the final sample of lenses. The dashed histogram shows the source redshift inferred from 
the COSMOS distribution. The thick histogram shows the final z s . p df, where the inferred source rcdshifts were weighted by the Sheth 
ct al. fitting function (multiplied by ct 8 ), as discussed in the text. For lenses with a measured z s , the grey vertical line shows the measured 
source redshift. Note that the PDFs have been normalized to unit area and that y-axis on the bottom row is on a different scale. 

those inferred by lensing and dynamical models, whereas 
the Chabrier IMF gives (on average) an underestimate. 
However, the IMF normalization may be mass depen- 
dent being more similar to a Chabrier IMF for the lower 
mass systems in the SLACS sample, with a <~ 200 km s _1 
(Barnabe et al. 2010; Auger et al. 2010; see also Cap- 
pellari et al. 2006, 2009). For simplicity, we adopt the 
Salpeter IMF normalization, although they can easily be 
converted to other normalizations, by multiplying by an 
appropriate factor (e.g. Auger et al. 2009). Total stellar 
masses for a Salpeter IMF are given in Table 3. 

LSD stellar masses were also calculated using the 
method described above with photometry from Koop- 
mans & Trcu (2002, 2003) and Treu & Koopmans (2004). 
This method could only be applied to 4 of the 5 LSD 
lenses, as photometry in multiple bands is required. Us- 
ing this method we inferred total stellar masses (in units 
of 10 n M©) of: 5.01±2.06, 2.24±0.58, 4.26±0.95 and 
2.49±0.82 for H1417+526, H1543+535, MG2016+112 
and 0047-281, respectively. 



5. LENS MODELING 

The modeling of our high resolution HST images is 
described in detail in Paper I. Here, we briefly outline 
the procedure followed and describe the relevant output 
parameters that we use in our joint analysis. Following 
standard image reduction, the bright lens galaxies were 
subtracted from postage stamp images of the lenses, iter- 
atively masking the lensed arcs. The remaining residual 
image was fitted using a flexibly parametrized source (an 
elliptical exponential profile with free position, ellipticity, 
orientation, size and flux) traced through a simple ellipti- 



cally symmetric lens potential. We model each lens mass 
distribution using a Singular Isothermal Ellipsoid (SIE) 
model (Kormann et al. 1994), with centroid fixed at the 
centroid of the galf it model lens light distribution. The 
lens has three free parameters: Einstein radius i?Ein, ori- 
entation cf> and axis ratio q. In all cases the lens redshift 
is known spectroscopically, allowing us to work with Ein- 
stein radius in kpc. 

The parameters of each lens were inferred by exploring 
their posterior PDF using a Markov Chain Monte Carlo 
(MCMC) sampler. While the samples do characterize 
the posterior PDF as required, the simplicity of the lens 
model and the idealized source light distribution mean 
that, in practice, the posterior width does not provide 
a good estimate of the uncertainty on, for example, the 
Einstein radius. Instead, we use independent reconstruc- 
tions in several bands (which we have for a subset of the 
lenses) to estimate the uncertainties (see e.g. Marshall 
et al. 2007), which are found to be about 5% on i?Ein- 

The properties of the lens models are discussed in de- 
tail in Paper I; here we are concerned only with the 
combination of the lensing mass with the dynamical 
mass from our spectroscopic measurements. With this in 
mind, the definition of elliptical radius in the lens model- 
ing was chosen to preserve the enclosed mass with circu- 
lar apertures, adopting the same definition as the SLACS 
survey (Koopmans et al. 2006; Bolton ct al. 2008b; Auger 
et al. 2009). The mass enclosed within the Einstein ra- 
dius, Msin, is essentially independent of the profile of 
the lens density profile: although we used an SIE model 
to infer Meih, we can treat it as a profile-independent, 
circularly symmetric enclosed mass and then combine it 
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Fig. 6.— Redshift distributions of the SL2S, SLACS and LSD 
surveys. All three distributions have been normalized to one. 

with a dynamical spherical mass estimate in order to in- 
fer the total density profile slope (see Section 7.4 below). 
For more discussion of the accuracy of this assumption 
see Auger et al. (2009). 

Three systems J140614+520252, J220629+005728 and 
J221606-175131, were observed with LRIS before at- 
tempting to perform a lens model. But subsequently we 
were not able to find a satisfactory lens model and we 
failed in measuring their Einstein radius essentially be- 
cause HST imaging is not available for them yet. Deflec- 
tor photometry along with effective radius and redshift 
are nevertheless reported in Table 2. A final conclusion 
on their lensing nature is left for future work (see also 
Paper I for further discussion). 

6. CHARACTERIZATION OF THE SL2S SAMPLE 

In this section we explore the basic properties of the 
SL2S lens galaxies, and compare them with two previ- 
ously studied samples of lenses. The SLACS sample is a 
lower redshift sample, 0.08 < z d < 0.51, of 85 lenses from 
the Sloan Digital Sky Survey (Bolton et al. 2006, 2008a). 
The LSD survey measured the internal kinematics of 5 
early-type lens galaxies over a large range of redshifts, 
0.48 < z d < 1.00, and masses (Treu & Koopmans 2004). 

6.1. The redshift distributions of SL2S lenses 

It is instructive to compare the distribution of lens red- 
shifts in the SL2S sample, described in the preceding sec- 
tion, with those of the lower and higher redshift reference 
samples, SLACS and LSD, respectively. In Figure 6, we 
show histograms of deflector redshifts for the SL2S (red), 
SLACS (grey) and LSD (black) samples. We see that, as 
anticipated, the SL2S lenses lie at higher redshift than 
the SLACS sample. The median and 68% confidence 
intervals are z d =0.494±g;^, 2 s =1.199±^g for the SL2S 
sample, z d =0.189t^o66 ^=0.606^;^7 for SLACS and 
3 d =0.810±g;g{$ z s =3.263i[^i for LSD. 

6.2. Comparison to previous lens samples 

The measured stellar velocity dispersions and effective 
radii provide us with a means to compare the SL2S, 
SLACS and LSD samples, to assess whether we are 
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Fig. 7. — Stellar velocity dispersions are plotted against the 
effective radius in the V band. The SL2S, SLACS (Auger et al. 
2009) and LSD (Treu & Koopmans 2004) data points are shown 
in red, grey and black, respectively. Note that R ef j is measured at 
the intermediate axis and the stellar velocity dispersions have been 
corrected to cr e 2- 
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Fig. 8. — The effective radius is plotted against the Einstein 
radius. The SL2S, SLACS and LSD points are shown in red, grey 
and black, respectively. Typical uncertainties on i?Ein an d R a ff are 
5% and 10%, respectively. 

studying the same types of galaxies. In Figure 7, we plot 
our spectroscopic velocity dispersions, a against the ef- 
fective radii, i? e ff inferred from the CFHTLS data above, 
in order to compare the types of galaxies selected in the 
SL2S, SLACS and LSD samples. The stellar velocity dis- 
persions are consistent between the samples, while the 
effective radii of the SLACS galaxies are slightly higher. 
Note that the stellar velocity dispersions have been nor- 
malized to a standard aperture, a C 2 for all samples. 

To compare the samples quantitatively, we applied the 
2D Kolmogorov-Smirnov test in the a — R c g plane. We 
found the KS statistic to be 0.28, and the significance 
level to be 91%. Given this high significance, we con- 
clude that we are drawing from the same population of 
galaxies, despite the marginal tendency towards larger 
sizes in the SLACS sample. 
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Figure 8 shows the distribution of R c g and i?Ein for 
the SL2S, SLACS and LSD samples. The higher red- 
shift samples (SL2S and LSD), have much higher physical 
-REin values while the effective radii are slightly smaller, 
as discussed above. The difference in i?Ein is partly due 
to the method used to select the lens candidates: re- 
solvable SL2S ring radii of 0.5-2 arcsec correspond to 
larger physical radii than those of the SLACS sample. 
SLACS lenses were pre-selected to be bright SDSS tar- 
get galaxies. This selection method ensures that they 
lie at relatively low redshift, which in turn means that 
the characterizing angular size of Einstein rings trans- 
lates to a smaller physical Einstein radii for the SLACS 
lenses than the SL2S lenses. Furthermore, the source 
redshifts are on average significantly higher for the LSD 
and SL2S samples, implying larger Einstein Radii for the 
same mass distribution. 

We conclude that we do seem to be measuring similar 
types of galaxies, with similar sizes, in a similar mass 
range: the SLACS and SL2S lenses are consistent with 
having been drawn from the same population. However, 
the different lens geometry (due to the difference in red- 
shift between the samples) means that we are sampling 
the mass density profile at a larger radius, which is ap- 
proximately a factor of 2 times greater compared to the 
SLACS lenses. 

7. THE INTERNAL MASS STRUCTURE OF SL2S LENSING 
GALAXIES 

In this section we use our new spectroscopic measure- 
ments to carry out a joint lensing and dynamics analysis 
of the SL2S sample, estimating the total mass density 
profile slope and dark matter fraction of each individual 
lens galaxy. We first describe the simple model within 
which we work, and then present our inferences of these 
two key parameters. 

7.1. Power-law density profiles 

We choose to work in the context of a spherically- 
symmetric power law total mass density profile model 
(e.g. Treu & Koopmans 2002; Koopmans et al. 2006; 
Suyu et al. 2010): 
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As in the lens modeling, we use capitalized R to denote 
projected radius. We can normalize these power law pro- 
files in terms of M^in, the robustly-estimated, profile- 
independent Einstein Mass, from Section 5 above: our 
approach is to use the well-measured Einstein radii and 
Einstein masses from the lens models of Paper I, and re- 
interpret them within the context of this spherical power- 
law model. 

The spherical power law model allows us to predict the 
observed spectroscopic velocity dispersion of the tracer 
stellar population, via the Jeans equation (see, e.g. Treu 
& Koopmans 2004; Koopmans et al. 2006; Auger et al. 
2010; Suyu et al. 2010). Normalizing the surface density 
and integrating to the (observable) Einstein radius i?Ein, 
we find that 
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Since M E i n = 7i"-R|in^crit, we can choose the two pa- 
rameters of our model to be Msm and 7'. In the next 
subsection we outline how the stellar velocity dispersion 
is predicted from these parameters. 

7.2. Stellar dynamics analysis 

As described by Suyu et al. (2010) and Auger et al. 
(2010), we predict the stellar velocity dispersion at the 
appropriate aperture radius using the spherical Jeans 
equation. Qualitatively we use the power law profile to 
compute the total spherical mass enclosed within radius 
r, and assume the tracer stars follow a Hernquist dis- 
tribution with scale radius related to the measured ef- 
fective radius, and then integrate to predict the stellar 
velocity dispersion, cr(M E i n ,7'). We direct the reader to 
Section 2.3 of Suyu et al. (2010) for all relevant equa- 
tions. We assume isotropic orbits for all of the systems, 
which is approximately found to be the case from a more 
detailed analysis of the resolved kinematics of several 
SLACS lenses (Barnabe et al. 2009). 

7.3. Joint analysis 

We propagate the uncertainties in the Einstein mass 
(which can be significant, for systems with no spectro- 
scopic source redshift), and the spectroscopic velocity 
dispersion, in the following way. The output of the lens 
modeling procedure for each lens is the posterior proba- 
bility distribution for the Einstein mass given the HST 
image data d, Pr(MEi n |d, z<j, z s ), characterized as a set of 
sample ME; n values - this can be viewed as the prior PD 
for Msin for the joint analysis. The analysis of that lens' 
spectrum yields the likelihood, Pr(cr obs |MEi n , 7'): which 
we assume to be Gaussian in cr obs , with mean equal to 
the Jeans prediction a (MEin, Y) from the previous sub- 
section. 

The posterior PDF that we seek for each lens is then 
Pr(M Ein , 7V bs , d, z d , z s ) oc Pr(a obs |M Ein , 7') 

•Pr(M Ein |d,z d ,z s ) 
■Pr(Y). (10) 
Since we are only working in the context of a single 
model, we do not compute the evidence to normalize the 
right-hand side of this equation; all we need are samples 
drawn from the posterior PDF Pr(M Ein , y\a° hs ). 

For the prior on the slope 7' we assume a uniform dis- 
tribution with limits of —1.2 and —2.8, which we find to 
be broad enough to enclose all the likelihood (and are in 
fact close to the mathematical limits required for the nor- 
malizability of the profile) . We draw an equal number of 
sample 7' values from this uniform distribution to match 
the AfEin samples, and then use the likelihood evaluated 
at each 2-dimcnsional sample position as weights: for 
each sample {MEi n ,7'} we solve the Jeans equation to 
calculate a(MEin,j'), and then evaluate the Gaussian 
function Pr(cr obs |M E i„, 7')- (For more details on impor- 
tance sampling, see the appendix of Suyu et al. 2010.) 
This set of weighted samples can then be used to com- 
pute integrals over the posterior, confidence limits, his- 
tograms to represent the marginalized distributions and 
so on. We present these inferences in the following two 
sections, and summarize our numerical results in Table 3. 
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Fig. 9. — Posterior probability distributions for 7' using a uni- 
form prior. The solid and dashed curves show 7' distributions 
for individual lenses with and without measured source redshifts, 
respectively. The shaded region indicates the posterior PDF for 
the mean of the Gaussian distribution from which the sample was 



inferred to have been drawn: {7') = 2.16 
tributions have been normalized to unit area. 
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Fig. 10. — The slope of the density profile is plotted against 
Einstein radius, projected mass within the Einstein radius, stellar 
cllipticity and effective radius. The error bars show the 16 th and 
84 th percentiles. The solid line shows the linear best fit to the data 
and the dashed lines indicate the scatter. 

7.4. The density profile slope 7' from lensing and stellar 

dynamics 

In total we have 11 SL2S lenses with both lensing and 
dynamical mass estimates. In Figure 9 we show the pos- 
terior probability distributions for the logarithmic den- 
sity profile slope 7', resulting from the joint lensing and 
dynamics analysis. We see that the mean of the popu- 
lation lies at (7') = 2.161q'q9, an< ^ * s shown by a shaded 
region in Figure 9. The intrinsic (Gaussian) scatter of 
the sample, inferred assuming Gaussian errors on the in- 
dividual 7' values, is Sy =0.25l ; 07 . For each of the 11 
lenses in the final sample, the median 7' (with 16 th and 
84 th percentiles) is given in Table 3 

We explore how 7' varies with: Reas-I RcS, total mass 



within i?Ein, axis ratio and effective radius of the lenses 
in Figure 10. In each panel, the solid line shows the 
linear best fit to the data and the dashed lines show the 
scatter. The gradients of the fits, including 16 th and 84 th 
percentiles are: 

di/d{R mn /R eS ) = -0.01±g;?i 

dj'/dM Retc = -0.03t - 17 



dj' /dq* - 
dj'/ dR e s ■■ 



-0.05 
n 22 +0.29 

■ u - lz -0.07' 



with Mjt eB in units of 10 M@, and R c s in kpc. The 
large uncertainties on the gradients indicate negligible 
trends between 7' and these variables. 

As discussed in Section 9, cosmic evolution of 7' could 
be mimicked by a dependence on i?Ein/-Reff- However, a 
negligible correlation between 7' and i?Ein/ R e s suggests 
that this is unlikely. 

We find no correlation between 7' and either R e s or 
the projected mass inside i? c ff, indicating that 7' is in- 
dependent of galaxy size and mass. 

7.5. The dark matter fraction in the SL2S lenses 

We now turn to our second key parameter: the dark 
matter mass fraction, /dm in the cores of lens galaxies. 
Following the SLACS analyses of Koopmans et al. (2006) 
and Auger et al. (2009), we use R e s/2 as the fiducial 
aperture radius for dark matter fraction estimation. The 
stellar masses calculated in Section 4.2 are total stellar 
masses. These must be corrected to the mass within our 
fiducial aperture by integrating (under the assumption 
that stellar surface density follows surface brightness) the 
de Vaucouleurs profile, which is given by: 



I(r) = /off e 



-7.67 



[te) 17 



(ii) 



The fraction of mass enclosed within R c g/2 is therefore 
given by: 



M.(iW2) _ Io R ° ie/2 I(r)rdr 
Mt : totai J °° I(r) r dr 



= 0.320. 



(12) 



We obtained total masses within the same radius by in- 
tegrating the power law total surface density profile: 

3-7' 

(13) 



M(R cS /2) fR eS /2 



Reii 



Total masses within the Einstein radius and total stellar 
masses are given in Table 3. The dark matter fraction is 
then defined by: 



/dm — 1 



M^Cgeff/g) 
M{R eS /2) ' 



(14) 



We do not just have a single number for each of 
M st (i? c ff/2) and A/(i? e ff/2); rather, we have probabil- 
ity distributions for each. We combine these to form 
the posterior PDF for /dm, by applying the formulae in 
Equations 12, 13 and 14 to each sample {M*,M} drawn 
from the product of Pr(M„) (from the SED modeling) 
and Pr(M ) (from the lensing plus dynamics joint anal- 
ysis) - since these PDFs are independent, the members 
of each {M*, M} pair can be drawn randomly from each 
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Fig. 11. — Posterior probability distributions for the dark matter 
fraction in the SL2S lenses. The dark matter distributions of lenses 
with measured source redshifts are shown in a solid line. The 
dashed lines show the distributions for lenses with no measured 
source redshift. Note that the distributions have been normalized 
to unit area. The shaded region indicates the posterior PDF for 
the mean of the Gaussian distribution from which the sample was 
inferred to have been drawn: (/dm) = 0. 42^ 



,+0.08 
-0.08- 



individual ensemble. The resulting PDF for /dm (visu- 
alized as a histogram of /dm samples) can have signifi- 
cant probability at /dm < 0, due to the uncertainty in 
each mass estimate, and the fact that our model only 
enforces positivity for the total and stellar mass distri- 
butions, not their difference. For consistency with the 
SLACS analysis, we allow negative values of /dm (equiv- 
alent to M* > M). We found that truncating the PDFs 
at /dm > for the SL2S lenses did not significantly affect 
the numerical results. 

In Figure 11 we show the resulting dark matter frac- 
tion posterior probability distribution for each SL2S lens, 
assuming a Salpeter IMF. The median and 16 th and 84 th 
percentiles of each dark matter fraction distribution are 
given in Table 3. We see that the mean of the popu- 
lation lies at 0.42^Q Qg, as shown by the shaded region. 

This distribution has an intrinsic width 0.201q'o7. Our 
results are sensitive to the choice of a universal Salpeter 
IMF in the stellar populations analysis. If we instead 
assert a universal Chabrier IMF, we find that the mean 
dark matter fraction is 0.68^Q Qg, while the scatter is es- 
sentially unchanged (0.11 tg^). As expected, a Chabrier 
IMF predicts a higher dark matter fraction, as shown in 
Table 3. 

In Figure 12 we show how the SL2S lens dark matter 
fractions vary with: i?Ein/^?cff, the total mass within 
i? c ff, axis ratio, density profile slope, and i? c ff- In each 
panel, the solid line shows the linear best fit to the data 
and the dashed lines show the scatter. The gradients of 
the fits, including 16 th and 84 th percentiles are: 
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Fig. 12. — Dark matter fraction plotted against REin/ Reff> t ne 
total mass within i? e fj, axis ratio, total mass density profile slope, 
and R c ff. The error bars show the 16 th and 84 th percentiles. The 
solid line shows the linear best fit to the data and the dashed lines 
indicate the scatter. 

with Mji c!t in units of \Q 11 Mq, R e g in kpc and /dm is the 
dark matter within R c g/2. There is evidence for a cor- 
relation between /dm and 7'. For all other variables, the 
large uncertainties on the gradients indicate that there 
is only a negligible trend with /dm- 

8. COSMIC EVOLUTION 

We now investigate evolution in the properties of mas- 
sive galaxies using our new lens sample. We focus on 
the two key quantities studied in Section 7: the density 
profile slope 7', and the dark matter fraction /dm- 

The low redshift reference measurements of 7' and /dm 
come from the SLACS analysis, which found density pro- 
files very close to isothermal inside one effective radius 
in a sample of 63 SLACS strong-lens early-type galaxies: 
(7slacs) = 2 - 078 ± °- 027 ' with a scatter of 0.16 (Auger 
et al. 2010). Likewise Auger et al. (2009) inferred a mean 
/dm within i? c ff/2 for 85 SLACS lenses of 0.3 with a scat- 
ter of 0.2, using a Salpeter IMF. These results are broadly 
consistent with our findings for SL2S systems. We can 
therefore anticipate only a mild cosmic evolution of these 
quantities. 

In Figure 13, we show how the density profile slope 7' 
varies with redshift, using the SL2S, SLACS and LSD 
samples together to cover the redshift range 0.05 to 1. 
We quantify this statement by fitting the 7'(.Zd) data 
with a linear relation in the mean slope, still including 
Gaussian scatter about that relation: 



(7')(^) = <7o) 



d(i) 



(15) 



For the SL2S 

d(i)/dz A 



data alone, we 

for the gradient and 



find ( 7o } = 2.22 +al7 



lfi+ - 48 

u - lo -0.51 



0.21! 
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Fig. 13. — Cosmic evolution of total mass density slope, 7'. 
The SLACS and LSD values were taken from: (Auger et al. 2010) 
and (Treu & Koopmans 2002; Koopmans & Treu 2003; Treu & 
Koopmans 2004), respectively. The error bars show the 16 and 
84 th percentiles. The best fit to the data is shown by the solid line 
and the scatter is shown by the dashed lines. 
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cr, (km s _1 ) 

Fig. 14. — The dark matter fraction within R c f{/2 as a function 
of measured velocity dispersion. The error bars show the 16 th 
and 84 th percentiles. For all samples, the velocity dispersion was 
normalized to a standard aperture, R e g/2. 

evolving 7' case, the scatter is Sy =0.23^oo6- When we 
include the SLACS and LSD data points, we find (jq) = 
2.12l° it d(y)/dz d = -0.25± f 2 , and S Y =0.17t° f 2 . 

These results are inconsistent with no evolution in the 
total density profile slope of massive lens galaxies since 
z ~ 1: the probability of the linear gradient in (7') being 
positive is just 2%. The lens data suggest (at approxi- 
mately the 2-cr level) that the mean total density profile 
of massive galaxies has become slightly steeper over cos- 
mic time. 

Figure 14 shows the dark matter fraction as a function 
of the measured stellar velocity dispersion, a e 2 for the 
SL2S, SLACS and LSD samples. Together, these samples 
cover a large rage in stellar velocity dispersion, however, 
no significant trend with central dark matter fraction was 



Fig. 15. — Evolution of dark matter fraction with redshift for a 
Salpeter IMF. The error bars show the 16 th and 84 th percentiles. 
The solid and dashed lines show the best fit to the data and the 
scatter, respectively. 

found. In Figure 15, we show how the dark matter frac- 
tion varies with redshift, using the SL2S, SLACS and 
LSD samples together to cover the redshift range 0.05 
to 1.0, and focusing on the Salpeter IMF. (The Chabrier 
IMF assignment only affects the overall normalization of 
the stellar masses and not their evolution.) Only 4 of the 
5 LSD lenses were included, since photometry in multi- 
ple bands was required to infer stellar masses (given in 
Section 4.2). We find that the mean dark matter frac- 
tion has not evolved strongly with cosmic time, however, 
there is marginal evidence for some change in the popu- 
lation. Again, we quantify this statement by fitting the 
/dm(^ci) data with a linear relation in the mean slope, 
still including Gaussian scatter about that relation: 



(/dm)^) = (/dM,o) H f DM ^ Zd ± Sf m 

OZ d 



(16) 



For the SL2S data alone, we find (/dm,o) = 0.381°,;^, 



d(f U M)/dz d = 0.07 



+0.36 
-0.36 



for the gradient and S/ D1 



0.18l£$ for the scatter. Including the SLACS and LSD 
data points, we find (/dm,o) = 0.27±°,;°J, d(f D M)/dz d = 
0.36±g;££, and S> DM = 0.13±g;g|. The probability of the 
linear gradient in (/dm) being positive is 98%: the 
lens data suggest (again at approximately the 2-cr level) 
that the mean projected dark matter fraction in massive 
galaxies, within half their effective radius, has decreased 
slightly over cosmic time. 

9. DISCUSSION 

By combining the inferred total mass density slopes of 
the SL2S, SLACS and LSD samples, we found a tanta- 
lizing suggestion that 7' has become slightly steeper over 
cosmic time. 

Before interpreting this result, and its implications 
for our understanding of the formation and evolution of 
early-type galaxies, it is important to address a potential 
source of systematic error. As shown in Figure 16, the 
ratio of 7?Ein/-Reff increases with redshift, mainly because 
the physical size of i?Ein increases with redshift (as dis- 
cussed in Section 6.2), but also because of the increasing 
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Fig. 16. — Evolution of i?.Ein/-f?eff with redshift. The ratio of 
^Ein/ficff tends to increase with redshift. 

source redshift. Since the slope of the total mass density 
is determined by measuring the total mass at two differ- 
ent radii, one of which is the Einstein radius, SL2S lenses 
sample 7' at larger radii than SLACS lenses. Therefore, 
a trend in i?Ein/ RcS with redshift could mimic the in- 
ferred evolution of 7'. However, for this effect to mimic 
the evolution of 7', the density profile would have to be- 
come shallower and then steeper with increasing radius, 
which is unlikely given the total mass density profile of 
local early-type galaxies. We also note that 7' and /dm 
vary differently with cosmic time, which is unexpected 
if the evolution being mimicked by a trend in R^m/RcS 
with redshift. As discussed in Sections 7.4 and 7.5, Fig- 
ures 10 and 12 show that both 7' and /dm show a negli- 
gible trend with i?Ein/-Reff, suggesting that the evolution 
of 7' is not an artifact due to a dependence on i?Ein/-Reff- 

For these reasons we conclude that true evolution of the 
average mass density profile is the most likely interpreta- 
tion of our findings, although additional measurements 
of 7' at different radii are needed to conclusively rule out 
alternate interpretations. 

It is important to stress that both the mean value 
and evolution of 7' are completely consistent with previ- 
ous results, (7o,slacs+lsd) =2.10±0.07and d(j')/dz d = 
—0.23 ±0.16 (Koopmans et al. 2006). It is only by virtue 
of our larger sample that we have been able to reduce 
the error bars and find marginal evidence for evolution. 
If this evidence for evolution is confirmed by larger sam- 
ples, it would indicate that growth of massive galax- 
ies since z ~ 1 has not occurred through dry (dissipa- 
tionless) mergers alone, since they preserve 7' (see e.g. 
Nipoti et al. 2009b). In dissipative merging events, 7' 
increases as a result of baryons cooling and sinking to- 
wards the center. Therefore, dissipative processes would 
be required to contribute, at least partially, to the evo- 
lution of early- type galaxies since z ~ 1. Although they 
cannot be the dominant process, because of tight limits 
on recent star formation in massive early-type galaxies 
since z ~ 1 (e.g., Treu et al. 2005), our result seems to 
suggest that they cannot be completely neglected either. 
This is consistent with evidence for a "frosting" of recent 
stars found in detailed studies of the stellar populations 



of massive galaxies (e.g., Trager et al. 2000). 

Another notable feature of the SL2S 7' distribution is 
its large intrinsic scatter. The intrinsic scatter in the 
SLACS sample is just 0.16±0.02 (Auger et al. 2010), 
compared to 0.25t";„" for SL2S. From the LSD lenses, 
Treu & Koopmans (2004) also found that the intrin- 
sic scatter in 7' was larger at z = 1. This may in- 
dicate an overall trend toward more complete dynam- 
ical relaxation over Gyr timescales, and perhaps a re- 
duced contribution of external convergence, due to line 
of sight structure in nearby systems, where the Einstein 
radii are smaller. Again, a larger sample size is re- 
quired to confirm this interpretation. Finally, we find 
that the dark matter fraction within half of the effec- 
tive radius has decreased slightly with cosmic time. This 
is again consistent with some contribution from dissipa- 
tional processes in early-type galaxy formation and evo- 
lution, where baryons move to the central region as they 
cool. However, it is important to keep in mind residual 
uncertainties in estimating the stellar mass, which may 
be redshift dependent, since the average age of the stellar 
populations is a function of cosmic time. Cosmic evolu- 
tion of the IMF would be another source of ambiguity in 
interpreting these results (van Dokkum 2008; Treu et al. 
2010). 

To conclude, we have developed a new method to 
estimate source redshifts, in the absence of spectro- 
scopic measurements. The method is completely gen- 
eral, does not require accurate multiband photometry of 
the sources, and allows us to infer 7' and /dm with errors 
that are just 2-3 times as large as if we had spectroscopic 
redshifts. Although spectroscopic redshifts are in gen- 
eral preferable, this may be a good strategy for exploit- 
ing future surveys of thousands of strong lenses, where 
wholesale spectroscopy of the complete sample may not 
be practical. 

10. SUMMARY AND CONCLUSIONS 

New spectroscopic measurements with deflector red- 
shifts and velocity dispersions were presented for 11 
lenses. These spectroscopic measurements were com- 
bined with lens models and photometry, described in 
Paper I, to infer the total density slope and dark mat- 
ter fraction of each of the 11 SL2S galaxies in the final 
sample. 

The main results are summarized below: 

1. The SL2S sample has a median deflector redshift, 
Zd = 0.494, source redshift, z s — 1.199 and velocity 
dispersion, <7 C 2 = 273 km s _1 . 

2. The SL2S, SLACS and LSD lenses are the same 
types of galaxies, however, the physical size of i?Ein 
is generally larger for higher redshift deflectors. 

3. We developed a new method to estimate the source 
redshift probability distribution function for lenses 
with no spectroscopic z s . This lack of accurate 
source redshift produces uncertainties on 7' and 
/dm that are only a factor 2-3 greater. Uncertain- 
ties would eventually further decrease with better 
multiband source photometry, and this would al- 
low large samples of high redshift lenses to be an- 
alyzed in the next generation of cosmological sur- 
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veys, where spectroscopy of all the systems may 
not be affordable. 

4. The average total density slope measured from the 
SL2S sample alone is: (-y' ) = 2.16^o.o9> w ^ n a scat- 
ter of 0.25+°;q°. 

5. Combining the SL2S 7' measurements with pre- 
vious analyses from SLACS and LSD, we find, 
(Y ) = 2.12±g;gt, d{i)/dz A = -0.25t°;i°, and S Y = 
0.17^q;q2. This suggests (at approximately the 2-cr 
level) that the mean total density profile of massive 
galaxies has become slightly steeper over cosmic 
time. 

6. Stellar masses were estimated using CFHT pho- 
tometry, and enabled us to disentangle the total 
and stellar masses. 

7. From this we inferred the dark matter frac- 
tion within half the effective radius: (/dm,o) = 
0.38^50, with a scatter S> DM = 0.18i£^ for the 
SL2S sample, using a Salpeter IMF. 

8. The combined dark matter fractions from all three 
samples suggest that the dark matter fraction 
within i?eff/2 has decreased slightly since z ~ 1. 
We find the mean dark matter fraction within 
i? e ff/2, gradient of evolution over cosmic time and 
scatter to be: (/ DM ,o) = 0.27±g;gg, d(fuu)/dz d = 
0.361°;^, and Syb M =0.13tjj;{g, respectively. 
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Note. — All spectroscopically-observed lens systems are shown. R e g values are measured at the intermediate axis. The signal to noise ratio per 
pixel was calculated over the rest wavelength range 4000-5000A. Typical uncertainties on i?Em and 7? e ff are 5% and 10%, respectively. The column, 
Flag, indicates whether HST or CFHT data was used to measure 7?Ein and R c g. 
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TABLE 3 

Inferred SL2S galaxy-scale lens properties 
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Note. — Only the 11 modeled lenses with measured velocity dispersions are shown. cr C 2 is the measured 
stellar velocity dispersion corrected to a standard aperture, asm values calculated from measured source 
redshifts are given without uncertainty and usie values calculated from z SjP df are given with the 16 th and 
84* percentiles (and were also weighted by the Sheth et al. fitting function and selection function). M* is 
the total stellar mass. M E i n is the total mass enclosed at the Einstein radius. Dark matter fractions are 
given at R c g/2 for both Salpeter and Chabrier IMFs. 



